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ABSTRACT 


Result!  of  >a  cotnj^frehensive  investigation  of-  Minuteman 
exhaust  plume  electrical  properties  are  presented.  Electron 
density  and  electron ?neutral  collision  frequency,  profiles  are  give^ 
for  the  Stag^  ^  plume  at  100  kft  and  the  Stage  2  plume  at  118,200 
and  30.p.'  kft.  This  work  is  motivated  by  the  need  to  determine  plume 
electrical  conductivities,  input  data  toan  evaluation  of  EMI?/exhaust 
'|rlume  interactions*  The  calculations  are  described,  in  detail  and 
include  analyses  of:  (i)  two>phase  noszle  flow  via  the  method  of 
characteristics,  (ii)  nonequilibriuni  chemistry  along  nozzle  and 
plume  streamlines,  (iii)  inviscid  plume  properties  via  the  method  of 
characteristics,  (iv)  turbulent  mixing  /  nonequilibrium  chemistry 
along  c'.i'"ved  plume  boundaries,  and  (v)  nonequilibrium'  chemistry 
behind  the  Mach  disc.  The  results  include:  gas  and  particle  prop¬ 
erties  and  chemical  kinetic  distributions  ofcall  neutral  and  charged 
species  concentrations  in  the  nozzle:  inviscid  boundaries,  shock 
structure,  Mach  number  contours  and  electron  densities  In  the 
plume* 
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NOMENCLATURiS 


A*  nozzle  throat  area 

Aex  nozzle  exit  area 

Cq  specific  heat  at  constant  pressure 

D  Mach  disc  diameter 

dj,^  nozzle  exit  diameter 

k,  ki  constants --defined  by  Eq.  (1), 

M  Mach  riutnber 

n  defined  by  p.  ^ 

n^  electron  density 

Pr  Pvandtl  number 

p  static  pressure 

r  radial  distance  from  centerline 

nozzle  exit  radius 

T  temperature 

W  mixture  molecular  weight 

X„-  electron  mole  fraction 

e 

X  axial  distance  from  nozzle  exit 

Y  mass  fraction 


Greek  Letters 

•y  ratio  of  specific  heats 

p  viscosity 

Vg  electron-neutral  collision  frequency 

Subscripts 

eq  equivalent  nozzle 

ex  nozzle  exit  plane 

g.  gas 

w  free  stream 

* 

t  liquid 

s  solid 
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I.  INTRODUCTION 


Interaction  of  an  electromagnetic/ pulse  (EMP),  following  a  nuclear 
explosion,  with  a  missile  in  flight  can  induce  skin  currents  which  seriously 
affect  the  missile's  electrical  system.  Because  the  exhaust  plume  is  itself  an 
electrical  conductor  currents  induced  in  the  plume  can  be  transferred  to  the 
skin  of  the  missile  via  its  bas:^.  Thus  an  important  problem  is  the  extent  of 
the  electrical  coupling  between  the  missile  body  and  exhaust  plunie.  In  order 
to  estimate  the  magnitude  of  this  problem  the  exhaust  plume  electrical  conduc¬ 
tivity  must  be  known. 

Minute  man  EMP/plume  interactions  are  under  investigation  at  the 
Air  Force  Weapons  Laboratory  (APV/L).  Their  approach  consists  of; 

(i)  design  of  an  electrical  circuit  which  (when  connected  to  the  missile)  will 
"simulate"  the  induced  skin  currents*  which  result  from  plume/body  electrical 
coupling,  (ii)  measurement  of  induced  skin  currents  and  (ii)  interpretation  of 
the  results  in  terms  of  the  Minuteman  exhaust  plume's  influence  on(ihe  skin 
currents  induced  by  EMP  throughout  the  flight  trajectory. 

AeroCi.em's  contribution  to  this  program  has  been  to  predict  Minute- 
man  exhaust  plume  electrical  properties  (electron  density  and  electron-neutral 
collision  frequency)  for  a  limited  number  of  first  and  second  stage  trajectory 
points;  these  predictions  have  been  used  (by  AFWL)  to  calculate  the  plume/ 
body  electrical  coupling.  Although  some  predictions  of  Minuteman  exhaust 
plume  electrical  properties  have  been  reported^  there  has  been  no  attempt,, 
prior  to  the  present  effort,  to  apply  the  available  systematic  proceduriis  (see, 
e.g.  Refs.  2-5)  for  calculating  plume  electrical  properties.  It  is  the  purpose 
of  this  report  to  present  the  results  of  our  comprehensive  investigation  of 
Minuteman  exhaust  plume  electron  densities  and  collision  frequencies --in 
which  we  utilize,  wherever  possible,  state-of-the-art  computational  techniques 
for  both  nozzle  and  plume  flows.  A  detailed  study  of  the  nozzle  flow  is  essen¬ 
tial  for  Minuteman  because  (i)  the  high  particle  loading  (»  30  >  by  weight) 
strongly  influences  the  gas  properties  and  composition  at  the  nozzle  exit  pUne 
and  (ii)  at  the  altitudes  considered  here, in  (100,  200,  300  kft)  the  only  source  of 
ionization  (upstream  of  the  plunie  shocks)  is  in  the  combustion  chamber  (i.e. 
there  is  no  afterburning);  thus  electron  densities  within  the  plume  are  largely 
controlled  by  nozzle  cheihical  kinetic  processes. 

The  systematic  calculation  of  Minuteman  plume  electrical  properties 
involves  a  series  of  individual  steps,  starting  with  the  following  data  (i)  propel¬ 
lant  composition  and  alkali  metal  impurity  levels,  (ii)  combustion  chamber 


*  This  equivalent  circuit  is  , known  as  a  plunoe  "simulator". 
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pressure,  (iii)  nuzzle  design  and  (iv)  fUght  trajectory.  These  steps  Include: 

1.  Thermochemical  equilibrium  calculation  of  comisustion  chamber  » 
properties. 

2.  Two -phase  method  of  characteristics  calculation  of  nozzle  gas  and 
particle  flow  properties. 

3.  Calculations  ofmeutral  and  charged  species  concentrations  along 
nozzle  gas  ‘  treamlines  via  a  nonequilibri 'm  chemistry  program. 

4.  Plume  method  of  characteristics  calculations  up  to  the  first  Mach 
disc  to  deterniine  the  plume  boundary,  inviscid  flow  properties 
and  'Shock  structure. 

5.  Neutral  and  charged  species  concentration  calculations  along  plume 
streamlines  via  a  nonequilibrium  chenustry  program. 

6i  Turbulent  mixing/nonequilibrium  chemistry  calculations  along 

curved  plume  boundaries  to  give  flow  properties  and  composition  in 
the  mixing  regions. 

7.  Preliminary  considerations  of  flow  behind  the  first  Mach  disc  to 
estimate  plume  electron  densities. 

Alternative  approaches  to  these  predictions  are  possible  although  few  are  as 
comprehensive  as  the  procedure /outlined  above.  These  are  mentioned  through¬ 
out  the  following  sections,  tpgetlier  with  our  reasons  for  their  rejection  in 
favor  of  the  present  approach. 

The  pertinent  results  for  (and  the  uncertainties  in)  each  phase  of  the 
calculations  are  discussed  below.  These  results  will  serve  as  a  good  start¬ 
ing  point  for  more  sophisticated  models  and  will  provide  information  for  those 
interested  in  Minuteman  nozzle  and  exhaust  plume  flows  for  other  applications. 
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II.  COMBUSTION  CHAMBER/EQUILIBRIUM  NOZZLE  FLOW 


Minuteman  first  and  second  stage  propellant  compositions  and  nozzle 
parameters  pertinent  to  a  one -dimensional  equilibrium  calculation  of  nozzle 
flow  properties  and  composition  are  given  in  Table  1,  No  measurements  of 
alkali  metal  impurity  levels  in  actual  Minuteman  propellants  are  available. 
However  AeroChem  measured^  antounts  of  Na  and  K  in  samples  of  the  Thiokol 
TPH-lCll  propellant  used  in  200P-lb  thrust  motors.  Assuming  that  the  impu¬ 
rity  content  of  a  propellant  is  primarily  dependent  on  the  ingredients --rather 
than  on  the  manner  in  which  a  particular  rocket  motor  is  manufactured --we 
can  safely  apply  these  measured  values  to  the  Minuteman  Stage  1  propellant. 
Using  similar  reasoning,  and  noting  that  the  dominant  source  of  impurities  is 
introduced  via  ammibnium  perchlorate/  we  assume  the  impurity  levels  to  be 
the  same  for  the  Stage  1  and  Stage  2  propellants.  The  uncertainties  in  the 
results*  due  to' this  assumption  are  expected  to  be  small. 

Tables  II  and  III  give  the  one -dimensional  equilibrium  chamber  and 
(Uozzle  flow  properties  for  Stages  1  and  2,  respectively,  computed  via  a  stan¬ 
dard' thermochemical  equilibrium  code.*  The  following  points  should  be  noted: 


(1)  The  quantities  of  combustibles,  CO  and  Hj,  in  the  exhaust 
products  are  very  high  (about  60  mole .%  for  Stage  I  and  about 
50  mole  %  for  Stage  2)  at  their  one -dimensional  nozzle  exit  pres¬ 
sures  (see  Table  IV).  There  is  therefore  a  large  potential  for 
exhaust  plume  afterburning. t 


(2)  Na  and  K  contribute  about  equally  to  the  total  positive  ion  concen- 
.tration  in  the  combustion  chamber.  Even  though  the  equilibrium' 
decay  of  (Na'*’!  is  larger  than  the  decay  of  [K‘*’J  in  the  nozzle  expan¬ 
sion,  both  Na^  and  will  be  included  in  the  nonequilibrium 
streamline  calculations  to  determine  their  relative  influence,  on 
electron-positive  ion  recombination. 


+  + 

(3)  The  Al-containing  positive  ions,  A1  and  AlOH  can  be  safely 

neglected  in  subsequent  kinetic  calculations  because  their  concen¬ 
trations  are  small  conq>ared  to  [K*];  also  [A102‘']  can  be  neglected 
because  it  is  snruU  compared  to  [Cl']. 


*  Throughout  this  report  the  phrase  ’'uncertainties  in  the  results"  refers  to 
uncertainties  in  exhaust  plume  electron  densities. 

i 

which. wiU  occiir  at  altitudes  less  than  about  75  kft 
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(4 )  About  30%  of  the  total  mass  of  the  mixture  is  in  the  form  of 
A^Qi  liquid  or  solid  parities.  Thus  differences  between  gas 
and  particle  temperatures  and  velocities  at  the  nozsle  exit  plane 
are  expected  to  be  large. 

Table  IV  gives  the  one -dimensional  equilibrium  (gas/particle  and 
thermochemical)  exit  plane  conditions  for  the  area  ratios  given.in  Table  1. 
These  arc  useful  reference  values  for  assessing  the  influence  of  paiticle  ther¬ 
mal  and  velocity  lags  on  nozzle  gas  prop  ties. 

Input  data.to  the  two-phase  method  of  characteristics  nozzle  program 
(see  Section  111)  requires  properties  of  the  "gas  only"  at<the  temperature  that 
exists  for  gas/particle  equilibrium.  Average  specific  heats  and  tnolecular 
weights  of  the  gas  are  given,  in  Table  V. 


III.  NOZZLE  FLOW 


Some  possible  approaches  to  computing  flow  properties  and  composi¬ 
tion  in  the  MinuUman  first  and  second  stage  nozzles  (nozzle  contours  are 
shown  in  Figs.  1  and  2)  and  their  limitations  are  listed  in  Appendix  A.  For 
the  present  study  we  chose  to  use  an  analysis  in  which  (i)  the  two-phase  two- 
dimensional  nature  of  the  flow  is  accounted  for  via  a  method  of  characteristics 
(MOC)  solution  and,  (ii)  the  results  of  the  MQC  program  are  input  to  a  non- 
equilibrium  calculation  of  neutral  and  charged  epecies  concentrations  along 
nozzle,  streamlines . 


A.  Gas  and  Particle  Properties 

The  two-dimensional,  two-jdiase,  conetant  Y  method  of  characteristics 
(MOC)  program’  used  in  this  study  is  a  proven  tool  for  making  nozsle  perfor¬ 
mance  calculations  and  is  in  wide-spread  use.  The  analysis' and  connputer 
code  are  discussed  in  detail  iii  Ref.  9.  The  cods  requires  input  data 
on  the  gas  properties  (chamber  conditions,  average  values  of  V ,  specific  heat, 
ntMlecular  weight,  etc.),  particle  properties  (j^rtlcle  mass  fraction,  and  size 
distribution),  ahd  nozzle  iise  and  contour. 


Specific  input  data  for  the  present  calc^atiods  are  given  in  Table  VI. 
The. average  gas  properties  Sye  dditermii^  at  the  gaz/^s^icle  equilibrium 
temperature  (see  Table  V).  values -W^r,  >p.aiil'»  bre  Veaeoni^le**  for 
calculations  on  Mintttemim-ty|^  piropellaat  JbeniiHM^tipn  products;  is  Ihe 

mass  fraction  rStie  in  the  ■ccs^blist'ltn  7^  pUfCteia  siMr^smbution 

is  based  priniarily  oo'a,prlvafln  Tido^';persoiuiel  who 

were  invclvsd'iii  ^/ad|ibrtits(i|pi^:d,s|pfj^^ 
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Minuteman  first  mot  /rs.  These  particle  sizes  are  consistent  with  other 

measurements  on  large  motors.**  We  assumed  the  particle  size  distribu¬ 
tion  for  Stage  2  to  be  identical  to  that  for  Stage  1'‘.  (We  will  demonstrate  later 
that  reasonable  uncertainties  in  the  particle  size  distribution  do  not  seriously 
influence  the^gas  properties.) 

The  first  step  in  the  calculation  procedure  is  to  establish  an  initial 
data  line  (M2  1)  for  the  MOC  solution.  For  a  gas/particle  flow  this  line  is 
downstream  of  the  throat  (see  Figs.  1  and  2).  Flow  prbjperties  in/the  tran¬ 
sonic  region  are  obtained  from  a  onerdimensional  sink  flow  solution^- -assum¬ 
ing  a  conical  inlet  section  which  joins  smoothly  to  the  throat  radius  of  curva¬ 
ture.  The  two -phase  MOC  solution  starts  at  the  end  of  the  transonic  region 
and  inarches  downslreani  in  the  *'usual"  manner.** 

Figures  3-8  show  the  calculated  distributions  of  pressure,  gas 
temperature  and  gas  velocity  along  the  wall  and  centerline.  The  relative  values 
of  wall  and  centerline  properties  are  typical**  for  bell-shaped  nozzles.  Just 
downstream  of  the  throat  the  <flow  along  the  wall  expands  more  rkpidly  than  the 
flow  along  the  centerline  which  has  not  yet  been  influenced^'by  the  expansion 
waves  emanating  from  the  wall.  However,  far  downstream  the  pressure  on 
the  centerline  drops  below  the  wall  pressure  as  the  how  along  the  centerline 
still  experiences  the  influence  of  the  rapid  expansion  occurring  further  up¬ 
stream  along  the  wall.  Shocks  which  may  form  in  the  nozzle  cannot  be  handled 
by  the  MOC  program.  Instead,  as  noted  in  the  Stage  1  results,  the  effect  of  a 
strong  isentropic  c6mpr>ezsion  is  experienced  (axial  distance  from  throat  " 

17  in.  )  by  all  fluid  properties.  In  the  Stage  2  nozzle  the  compression 
«  37  in.  from  the  throat)  is  weak  aiui  only  the  gas  velocity  (Fig.  8)  is  signifi¬ 
cantly  affected. 

Particle  velocities,  temperatures  and  mass  flow  rates  for  each  size 
particle  (2,  4,  6,  and  8p)  are  shown  in  Figs.  9-14;  limiting  particle  stream¬ 
lines  (LPS)^  are  shown  in  Figs.  1  and  2.  Particles  are  seen  (Figs.  11  and  12) 
to  remain  at  their  melting  temperature  (2318  K)  for  a  considerable  portion  of. 
the  nozzle  flow.  It  is  also  interesting  to  note  that  the  2p  particles  on  th  ) 
centerline  are  sufficiently  small  to  respoxd  to  the  Stage  1  nozzle  compression- - 
as  evidenced  by  the  decreaee  in  velocity  (Fig.  9)  and  increase  iii  temperature 
(Fig.  11)  near  the  point  where  the  conapression  wave  intersects  the  axis.  Parti¬ 
cle  inass  flow  rates  (|^r  unit  .area),  which  are  needed  input  data  to  the  calcula¬ 
tion  of  heterogeneous  electron/positive  ion  recombination  ratea  are  presented 
in  Figs.  13  and  14. 


t  LP8  is  define  as  the  beua4iag  partioU  stredmliae  for  a  given  particle 
siee,  ii.e.  .fi|  phirtic^  of  Ih*  given  sthe. are  contained  between  the  center 
line  aBd'|hii-'|j(^. 


In>  order  t'i  determine  the  inflaence  of  uncerteintiee  in  the  particle 
•ize  distribution  oti  the  gas  properties,  edditiooal  (Stage  2)  nosale  calculations 
were  made  assuming  uniform  particle  diameters  of  3p  and  8p.  This  range 
ot  particle  diameters  appears  to  be  reasonable**  for  rocket  motors  approxi¬ 
mately  equal  in  size  to  the  Minute  man  motors.  The  results  of  these  para¬ 
metric  calculations,  together  with  results  using  the  standard  (2>  4.  6.  Sp) 
particle  size  distributiiim  are  shown  in  Figs.  15-17r-where  it  can  be  noted  that 
even  serious  errors  iu  particle  size,  distributions  (within  Ihe  3  -  8  p  range)  do 
not  significantly  influence  the  gas  properties.*  This  is  a  very  important  con¬ 
clusion  since  there  is  usually  a  large  uncertainty  in  the  actual  size  distribution 
generated  by  highly  aluminized  solid  propellant  systems. 

Nozzle  exit  plane  gas  and  particle  properties  (input  to  the  exhaust 
plume  MOC  program  and  the  turbulent  ndxing/nonequilibrium  chemistry  plunne 
program)  are  given  in  Figs.  18-21.  Radial  distributions  of  gas  properties 
(Figs.  18  and  20)  are  as  expected.  Properties  are  seen  to  be  nearly  constant 
from  the  centerline  up  to  "  0.6,  amd  then  vary  rapidly  up  to  the  nozzle 
lip.  Of  particular  importance  are  the  large  ditferences  between  actual  and 
one -dimensional  equilibrium  exit  plane  properties.  Use  of  the  one -dimensional 
properties  in  plume  calculations  could  therefore  result  in  serious  errors. 

Exit  plane  distribution  of  particle  properties  appears  to  be  consistent  with  the 
gas  properties,  it  is  interesting  to  no^e  that,  for  Stage  1.  (Fig.  19a)  both  6p 
and  8p  particles  remain  at  their  melting  temperature  across  the  entire  exit 
plane.  Also  note  that  there  is  little  difference  bstweeu  gas  and  particle  flow 
angles  (Figs.  19b  and  21b). 

Figure  22  shows  that  uncertainties  in  particle  size  distributions 
(over  the  3p  and  8p  range)  do  not  have  a  serious  effect  on  Stage  2  exit  plane 
properties. 

The  above  results  are  then  used  as  follows: 


(1)  Gas  pressure,  temperature  and  velocity  distributions  along  the  wall 
and  centerline  are  input  data  to  the  calculation  of  nonequillbrlum 
neutral  and  charged  species  concentrations  in  the  nozzle. 


(2)  Nozzle  exit  pUah  gas  pr^ertUs  are  input  data>to  the  exhaust  plum 
MOC  and  turbulent  ihixiag/aoneqniUbrium  chemistry  progranui. 


(3)  Particle  prope^^les  are  used  to  calculate  heterogeneous  elsctron^on 
recombinatiesi  ratas.  (There  was,  unfoirtittaately*  insufficient  time 
to  perform  Aesa  ealculatiOM  for  Resent  study.  Errors  in 
exhaust  pluihe  eleotron  fhmifUs  iai^eiaced.by  neglscting  heteroge¬ 
neous  recomMillfeai  a^^  ia  iectlMi  m.  B.  All  necsssary 

datdibr  'Inliiir#  calcdlild'lihi  iire/j|lbiii>tii  this  rf^rt.) 


*  The  heyjarsitiililr  )#' 


IS  Iractten  (T^Me  y)^ 
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B.  Neutral  and  Charged  Specio  Concent r>tion» 

A  new  computer  code  had:  to  be  written  to  make  ndnequilibrium 
calculations  of  neutral  and  charged  species  concentrations  along  streamlines 
in  a  two-phase  gas/particle  flow.  This  code  takes  1/iput  data  from  the  noz/.le 
MOC  calculations  and  integrates  the  species  continuity  equations  along  stream¬ 
lines.  An  implicit  numerical  integration  procedure  was  used  siniilar  to  that 
incorporated  into  the  JANNAF  ODK  program.*’* 

Figures:  3-8  give  the  input  gas  properties  for  the  Stage  1  and  Stage  2 
n\  zzle  wall  and  centerline  nonequilibrium  calculations.  Initial  conditions 
were  specified  at  the  beginning  of  the  transonic  flow  regHA--as  determined 
from  the  two-phase  MO^  program;  these  are  given  in  Table  VU.  The  initial 
mole  fraction  of  Na^,  K  ,  Cl*,  and  e"  were  adjusted  slightly^(in  proportion  to 
their  absolute  values)  to  maintain  charge  balance  because  Ai  ,  AlOH^  and  AIO2 
were  neglected. 


The  reaction  mechanism  and  rate.. coefficients  for  the  nozzle  stream¬ 
line  calculations  are  given  in  Table  VIU.  The  rate  coefficients  are  those 
recommended  by  Jensen  and  Kurtius*’  and  recently  updated  by  Jensen  and 
Jones;  they  are  considered  to  be  the  best  available  values.  The  wide  range 
of  the  error  bounds  for  many  reactions  (particularly  those  directly  affecting 
charged  species  concentrations)  must  be  accounted  for  in  assessing^the 
accuracy  of  results  (this  point  will  be  discussed  later). 

Figures  23-30  give  species  mole  fraction  distributions  along  the  wall 
and  centerline  for  the  Stage  1  and  Stage  2  nozzles  .as  conoputed  by  the  new  non¬ 
equilibrium  streamline  code.  The  following  points  should  be  noted: 

(1)  The  major  neutral  species  mole  fractions  remain  approximately 
constant  (Figs.  23  ai^  27). 

(2)  H  and  Cl  mole  fractions  freeze  rapidly,  while  the  other  free  radicals 
(O,  OH)  show  a  continuous  decay  (Figs.  24,  25,  28  and  29). 

(3)  NaCl  aixl  KCl  moie  fractions  remain  approximately  constant  while  Na 
and  K  mole  fractions  freeze  rapidly  (Figs.  25  and  29). 

(4)  Electron  mole  fractibos  decay  rapidly  to  their  minimum  values 
(Figs..  26  and  30).  The  subsequent  increase  in  (e*]  occurs  because 
^e  three -body  recombination  rates  (Reactions  17  and  18)  become 
smpli.  cotapared  to  ihb  rbversb  of  Reaction  22,  i.  e.  H  4  Cl*  HCl  4 
e',  which  domlnales  oyir  .the  forward  reaction  at  low  temperatures. 

•  ■  '  11  ■  Jr 
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This  iricrease  in  [e  ]  has  important  implications  ^ith  regard  to  the 
plume  calculations  (Section  iV.-B).  Note  that  [Na  ]  recombines 
much  faster  than  [K^]i  and  Na-cohtaining  species  can  safely  be 
neglected  in  subsequent  kinetic  calculations  of  this  type.* 

The  differences  between  wall  and  centerline  values  reflect  the 
differences  in  pressure/temperature  histories.  Only  [O],  [Oj]  and  [H]  are 
seen  to  have  significantly  large  radial  variations  at  the  exit  plane. 

The  influence,  on  electron  mole  fractions,  of  varying  the  rate 
coefficients  for  Reactions  (19)  and  (22)  (Table  VIII)  is  shown  in  Fig.  31.  We 
chose  to  vary  these  particular  rate  coefficients  because:  (i)  previous 
studies^^’^^  showed  plume  electron  densities  to  be  quite  sensitive  to  the  rate 
of  Reaction  (19) --and  using  ki^/lO  in  afterburning  plume  calculations  gave 
better  comparisons  with  attenuation  data  than  using  k|9  and  (ii)  [e  ]  is  known 
to  be  sensitive  to  the  rate  of  Reaction  (22)- -which  has  a  large  lower  bound. ^ 
The  results  of  this  parametric  study  (Fig.  31)  show  that  varying  ki9  and  k22 
by  factors  of  10  and  100  respectively  (which  are  within  the  error  bounds  cited 
in  Tabic  Vlll),  has  relatively  little  effect  on  [e'J  at  the  nozzle  exit  plane.  If 
heterogeneous  electron-positive  ion  recombination  rates  were  included  the 
exit  plane  values  of  [e’J  would  decrease  by  perhaps  a  factor  of  2- -making  the 
values  given  in  Figs.  26  and  30  conservatively  large. 


INVISCID  PLUME  PROPERTIES  AND  SHOCK  STRUCTURE 


:4f: 


The  method  of  characteristics  (MOC)  is  the  standard  technique  for 
computing  inviscid  underexpanded  rocket  plume  properties.  When  coupled 
with  the  Rankine -Hugoniot  equations,  which  describe  property  variations 
across  shocks,  it  is  an  extremely  useful  method  for  determining  inviscid 
plume  boundaries,  internal  shock  structure  and  flow  properties --up  to  the 
hrst  Mach  disc,  t  The  computer  code  used  to  make  these  calculations^  was 
developed  at  Lockheed,  and  has  been  demonstrated*^  to  give  predictions  of 
plume  boundaries  and  shock  structures  whicli  are  in  good  agreement  with 
experimental  observations. 


However  the  influence  of  Na  on  the  initial  positive  ion  concentrations 
will  have  to  be  accounted  for. 

The  recommended  rate  coetticient  for  Reaction  (22)  shown  in  Table  VUI 
is  larjj^er,  by  two  orders  bf  inagnitude  than  the  value  originally  recom¬ 
mended  by  Jensen  ai)4viCursittS,^*  , 

Calculations  ’dowi|Str«Nsm..6l  t|^‘'fisbt  |ds|c^^  diic  are  discussed  in  Section 
•YU. 


*  ■■  ‘v  «■  ‘  \  •!.  * ^3* . 't  » .1  /  K  ,,  K’  .  i 
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In  the  present  calculations  we  assume  that  the  particles  have  a 
negligible  influence  on  gas  properties  in  the  p'  .le.  This  appearr  to  be  a 
reasonable  assumption  since  the  gat?  expand^'  very  rapi  Uy  downstream  of  the 
exit  plane,  resulting  in  decreased  particle  drag  and  gas/particle  heat  transfer 
rates.  The  errors  introduced  in  the  calculations  of  plume  electrical  properties 
by  imposing  this  assumption  are  negligible. 


A.  Location  of  First  Mach  Disc 

.  position  of  the  first  Mach  disc  in  underexpanded  free  jets  has 
been  studied  experimentally  by  several  investigators.^*"*^  Lewis  and  Carlson** 
developed  a  semi-empirical  relationship  for  the  axial  position  of  the  Mach 
disc  based  upon  measurements  on  jets  from  a  15*  half  angle  conical  nozzle  into 
stUl  air  They  showed  that 


where  the  empirical  constants,  and  kj  were  found  to  be  0.69  and  0.5, 
respectively. 


Additional  studies**"*^  have  confirmed  the  form  of  Eq.  (1)  and  extended 
the  correlation  to  account  for  the  effect  of  a  supersonic  free  stream;  thus  for 
Mo,  >  I, 


x/dM  =  k  M  Y 
'  N  ex  'ex 


(2^ 


Equation  (2)  correlc  .es  da  a  over  a  wide  range  of  pressure  ratios,  jet 
exit  Mach  numbers  and  free  stream  Mach  numbers  and  can  therefore  be  used  to 
predict  first  Mach  disc  locations  in  rocket  plumes  with  a  relatively  high  confi¬ 
dence  level.  * 

The  influence  of  solid  particles  on  the  position  of  the  Mach  disc  has 
also  been  experimentally  studied  by  Lewis  and  Carlson.**  They  measured 
x/d^  as  a  function  of  particle  mass  fraction  and  found  the  following  correction 
factor  (to  Eq.  (2)), 


The  influence  of  nozzle  exit  flow  angle  on  first  Mach  disc  locations  has  not 
yet  been  experimentally  determined.  Equation  (2)  covers  the  range  of  flow 
angles  from  0  to  li>  deg,  which  is  the  range  of  interest  for  the  Stage  1  and 
Stage  2  nozzles. 


t 
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=  [ 


1  +  0. 197  M^j^*"** 


(3) 


Therefore,  for  gas-particle  flows  with  a  supersonic  free  stream  the 
position  of  the  Mach  disc  can  be  located  from 


0.69M  Y 


ex 

f(Yp) 


0.5 

ex 


(4) 


Experimental  data  on  the  Mach  disc  diameter  ar  '  given  by  Crist, 
et  and  D'Attorre  and  Harshbarger,^  but  these  data  have  not  been  reduced 
to  empirical  relationships.  It  has  however  been  found  that  the  diameter  has 
the  same  pressure  ratio  dependence  as  the  axial  location,  i.e. 

Because  of  the  lack  of  more  definitive  data,  we  determined  D  by  using  the 
MOC  solution  up  to  the  value  of  x/dj^  determined  from  Eq.  (4),  The  Mach  disc 
diameter  was  then  taken  to  be  the  intercepting  shock  diameter  at  the  axial 
location  of  the  disc. 

B.  Equivalent  Single  Nozzle  Approximation 

Minuteman  first  stage  consists  of  a  four-nozzle  assembly  (see  Fig. 
32).  Thus  the  interaction  of  four  exhaust  piumes  should,  in  principle,  be 
taken  into  account  in  calculating  first  stage  electrical  properties.  However 
techniques  for  making  interacting  plunne  flow  field  calculations  have  not  yet 
been  developed.  We  therefore  utilized  an  equivalent  single  nozzle  approxima¬ 
tion  in  which  the  four  nozzles  shown  in  Fig.  32  were  replaced  by  a  single 
nozzle,  with  the  same  total  mass  flow  rate  as  the  four  nozzles  (i.  e.  (^  M)eq  ~ 

2  r^).  At  values  of  x/r^  large  compared  to  the  plume  intersection  distance, 
this  assumption  should  be  reasonably  good.  However  the  assumption  breaks 
down  in  the  vicinity  of  the  missile  base  and  a  detailed  analysis  of  the  inter¬ 
secting  plume  flow  field  would  be  needed  for  an  accurate  description  of  the 
electrical  properties  in  this  region. 
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C.  Results 

Figure  33  shows  the  plume  boundary,  intercepting  shock,  Mach  disc 
and  constant  Mach  number  contours*  for  Stage  1  at  100  kft  (free  stream 
Mach  number  =  5. 15)  as  computed  from  the  Lockheed  MOC  program^®  (using 
the  input  data  given  in  Table  IX).  Note  the  rapid  expansion  in  the  region 
contained  by  the  shock  and  the  plume  centerline --so  that  M  >  8  before  the 
first  Mach  disc.  Also  observe  that  the  Mach  numbers  are  relatively  constant 
in  the  region  between  the  shock  and  the  boundary.  Figure  34  shows  the  influ¬ 
ence  of  free  stream  Mach  number  on  the  size  of  the  plume  and  position  of  the 
Mach  disc. 


The  results  for  Stage  2  at  118  kft,  200  kft  and  300  kft  are  shown  in 
Figs.  35-37,  respectively.  At  118  kft  the  ratio  of  nozzle  exit  pressure  to  the 
pressure  behind  the  external  shock  (not  shown  in  Fig.  35)  is  only  about  2 --too 
low  for  a  Mach  disc  to  form.  Instead  we  predict  that  the  diamond  pattern  (i.e. 
a  negligibly  small  Mach  disc)  typical  of  plumes  with  exit  to  free  stream  pres¬ 
sure  ratios  slightly  greater  than  «inity  will  form.  Also  shown  in  Fig.  35  are 
the  locations  of  the  streamlines  along  which  kinetic  calculations  were  made 
(see  Section  V).  The  Stage  2,  300  kft  results  (Fig.  37)  show  that  the  boundary 
and  shock  approach  each  other  with  increasing  altitude  and  the  Mach  disc 
moves  far  downstream  and  is  located  at  x/rj^  =  770. 


V.  NONEQUILIBRIUM  CHEMISTRY  -  PLUME  STREAMLINES 


The  nonequilibrium  chemistry  calculations  along  plume  streamlines 
were  made  with  the  same  computer  code  used  to  determine  neutral  and  charged 
species  concentrations  in  the  Stage  1  and  Stage  2  nozzles.  Input  data  to  the 
code  were  pressure,  temperature  and  velocity  distributions  along  individual 
streamlines,  as  defined  by  the  MOC  calculations,  and  the  species  mole  frac¬ 
tions  at  the  nozzle  exit  plane,  given  in  Figs.  23-30.  The  streamlines t  were 
followed  from  the  exit  plane  to  the  shock- -then  downstream  of  the  shock  to  the 
plume  boundary.^  Calculations  were  made  along  a  sufficient  number  of  stream¬ 
lines  to  describe  electron  densities  throughout  the  inviscid  plume;  some  of 


*  All  plume  properties  e.  g.  pressure,  temperature,  etc.,  can  be  computed 
from  the  standard  isentropic  flow  relationships  using  the  values  of  V  , 
stagnation  presikure  and  stagnation  temperature  listed  in  Table  IX. 

t  Only  streamlines  that  intersect  the  shock  upstream  of  the  triple  point 
were  considered. 

t 

t  The  flow  behind  the  (relatively  weak)  intercepting  shock  is  essentially 
frozen. 
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these  results  are  presented  in  Figs.  38-40.  Figure  38  shows  the  electron 
density  (and  electron  mole  fraction)  variation  along  the  centerline  of  the 
Stage  1  plume- -up  to  the  Mach  disc.  It  is  interesting  to  note  that  Xg- 
increases  -  -because  the  reaction  H  +  Cl  "*•  HCl  +  e"  is  still  significantly  fast 
at  temperatures  slightly  less  than  the  nozzle  exit  plane  temperature.  How¬ 
ever  X^-  rapidly  approaches  a  "frozen"  value,  after  which  ng  decreases  in 
proportion  to  the  gas  density.  Figure  39  shows  the  influence  of  the  value  of 
kzz  used  in  the  kinetic  calculations  on  Xg-  in  the  Stage  2  plume.  Using  kn 
recommended  by  Jensen  and  Jones^^  results  in  an  order  of  magnitude  increase 
in  Xg-  before  the  first  shock  intersection  with  the  plume  axis  (see  Fig.  35). 

This  increase  in  Xg  -  was  unexpected--since  it  has  been  generally 
thought  that  Xg_  decreases  in  inviscid  plumes  due  to  the  three -body  recombi¬ 
nation  reactions  (17)  and  (18)  (Table  VIII).  However  the  thermodynamics  of 
Reaction  (22)  favors  the  backward  reaction  as  temperature  decreases.  Thus 
Xq-  is  controlled  by  a  two-body  reaction  which  produces  electrons  rather 
than  a  three-body  reaction  which  depletes  electrons.  This  situation  occurs 
even  if  k22/l00  (which  is  within  the  error  bounds  recommended  by  Jensen  and 
Jones)  is  used  in  the  calculations  (as  seen  in  Figs.  39  and  40). 

Whether  Xg-  actually  increases  in  an  inviscid  plume  has  not  been 
chucked  out  experinientally.  It  is,  of  course,  possible  that  the  mechanism 
listed  in  Table  VIII  is  not  complete,  and  that  some  competing  two-body  recom¬ 
bination  reactions  have  been  omitted.  Jensen^®  has  suggested  that  the  forma¬ 
tion  of  hydrated  ions  (e.g.  K*H20'^)  in  the  chamber  and  the  subsequent  recom¬ 
bination  reaction,  K’H20'^  +  e  Products  might  influence  Xg-  in  the  nozzle 
and  plume  since  the  recombination  rate  for  this  two -body  reaction  is  faster 
than  the  three-body  recombination  reactions  (17)  and  (18).  However  there 
was  insufficient  time  in  this  short  study  to  explore  the  above  (and  other)  possi¬ 
ble  additions  to  the  reaction  mechanism  of  Table  VIII.  We  therefore  chose 
(as  a  compromise)  kzz/lOO  for  the  "standard"  electron  density  calculations 
reported  in  Figs.  46,  48,  50,  and  52  ,  and  estimate  these  results  to  be  within 
about  a  factor  of  3  of  the  actual  values. 

The  results  shown  in  Fig.  40  for  kjj/lOO  give  ng  distribution  in  the 
im/iscid  plume  (up  to  the  intercepting  shock  and  Mach  disc)  for  the  Stage  2 
altitudes  under  investigation. 
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VI.  MIXING/NONEQUILIBRIUM  CHEMISTRY  ALONG  CURVED 

PLUME  BOUNDARIES 


The  discussion  to  this  point  has  demonstrated  how  electron  densities 
in  the  inviscid  plume  have  been  computed.  The  next  step  is  to  couple  these 
results  with  a  calculation  of  the  properties  in  the  mixing  region.  An  analysis 
which  completely  couples  the  mixing/nonequilibrium  chemistry  equations 
(along  the  curved  plume  boundary)  to  the  method  of  characteristics  solution  is 
beyond  the  scope  of  the  present  study. For  the  present  calculations  we  have 
utilized  an  approximate  method  in  which  the  computed  pressure  distribution 
along  the  inviscid  plume  boundary  (which  takes  into  account  the  external  plume 
shock  caused  by  the  expanded  plume  boundary)  is  input  to  the  mixing/non- 
equilibrium  chemistry  program.^'^  We  then  determine  the  position  of  the  divid¬ 
ing  streamline  T  from  the  program  and  assume  that  the  plume  properties  com¬ 
puted  along  the  dividing  streamline  are  the  same  as  those  along  the  inviscid 
plume  boundary.  The  radial  property  distributions  and  width  of  the  mixing 
region  are  then  taken  directly  from  the  computer  output.  This  type  of  solutioi. 
is  reasonable  only  until  the  computed  mixing  region  intersects  the  plume  axis. 
Downstream  of  this  region  the  radial  property  distribution  must  be  estimated. 
Further  downstream  however,  the  effects  of  the  curved  boundary  become 
small,  the  plume  pressure  becomes  nearly  constant  and  equal  to  the  ambient 
pressure  and  the  axisyrnmetric  mixing/nonequilibrium  chemistry  program 
should  give  reasonably  accurate  plume  electron  densities. 

Input  data  to  the  mixing/nonequilibrium  chemistry  calculations  include 
the  species  mole  fractions  at  the  nozzle  exit^  given  in  Figs.  23-30  and  the  pres¬ 
sure,  tempero.ture  and  velocity  data  given  in  Table  X.  Figure  41  shows  the 
plume  boundary  pressures  computed  from  the  MOC  program,  while  Fig.  42 
shows  the  relative  position  of  the  dividing  streamlines  and  plume  boundaries. § 


k 
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*  Such  an  analysis  and  computer  code  will  be  developed  under  a  subsequent 
AFRPL  contract.  In  our  initial  work  on  this  problem  however  we  have 
formulated  the  appropriate  equations  and  have  considered  various  numerical 
solutions. 

t  The  dividing  streamline  is  the  streamline  which  starts  at  the  nozzle  lip. 

t  We  assume  that  the  species  mole  fractions  remain  frozen  at  their  exit 

plane  values  during  the  Prandtl-Meyer  expansion  from  the  nozzle  lip  pres¬ 
sure  to  initial  boundary  pressure. 

§  No  results  are  shown  for  Stage  2,  300  kft  since,  mixing  has  little  influence 
on  electron  densities  at  this  altitude. 
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P'or  Stage  1  (100  kft)  and  Stage  2  (118  kft)  the  new  superposition  technique 
described  above  appears  to  be  reasonable;  however  for  Stage  2  (200  kft)  the 
plume  boundary  location  is  much  different  from  the  dividing  streamline  loca¬ 
tion,  resulting  in  a  greater  inaccuracy  in  predicted  electron  densities. 
Temperature  and  electron  mole  fraction  distributions  along  the  dividing 
streamlines  are  given  in  Fig.  43. 

Ihe  mixing  region  was  assumed  to  be  tui’bulent  for  all  calculations. 
The  eddy  viscosity  model  used  was  that  of  Donaldson  and  Gray,^°  modified  by 
Tufts  and  Smoof'*  to  account  for  non-zero  free  stream  velocities. 


VII.  FLOW  PROPERTIES  BEHIND  FIRST  MACH  DISC 


There  has  been  no  published  work  to  date  which  gives  a  good 
computational  scheme  for  determining  flow  properties  and  composition  behind 
the  first  Mach  disc  (within  the  second  shock  cell).  We  have  made  a  liniited 
attempt  to  develop  such  a  computational  procedure  by  investigating  the  mixing 
that  occurs  along  the  slip  line,  between  the  fluid  passing  through  the  Mach 
clisc(normai  shock)  and  the  fluid  passing  through  the  oblique  shock.  (Thus  each 
fluid  lias  a  different  total  pressure.)  There  was  insufficient  time  however 
to  make  complete  calculations  with  the  AeroChem  Mi.xing/Nonequi)  ibrium 
Chemistry  program.*  Therefore  electron  densities  along  the  centerline  be¬ 
hind  the  Mach  disc  were  determined  from  the  streamline  kinetic  program. 

These  calculations  were  made  for  the  Stage  1,  100  kft  plume.  The 
initial  pressure,  temperature  and  velocity  were  computed  from  normal  shock 
relations;  initial  species  mole  fractions  were  assumed  to  be  equal  to  their 
values  upstream  of  the  Mach  disc  (see  Table  XI).  Pressure,  temperature  and 
velocity  distributions  along  the  centerline  (up  to  the  start  of  the  second  shock 
cell- -at  *  50)  were  estimated  by  reference  to  other  MOC  calculations. t 

These  estimates  are  given  in  Fig.  44.  Also  shown  are  initial  calculations 
made  with  the  mixing/nonequilibrium  chemistry  program.  Note  the  relatively 
large  changes  ui  velocity  and  temperature  due  to  mixing,  indicating  the  impor¬ 
tance  of  viscous  effects  in  this  region. 


*  This  will  be  pursued  in  a  future  contract  with  AFRPL. 

t  with  lower  nozzle  exit  to  ambient  pressure  ratios,  which  gave  results 
behind  the  first  oblique  shock  intersection  with  the  axis,  and  up  to  the 
second  shock  intersection. 


14 


TP-281 


Figui-o  45  shows  the  very,  large,,predicted  centerline  electron  mole 
fractions  and  electron  densities  downstream  o£  the  Mach  disc.  We  expect  that 
the  effects  of,  mixing  along  the  slip' line  will  be  to  lower  peak  electron  densities 
and  increase  the  electron  decay  rates.. 


VIII.  ELECTRON  DENSITY  AND  COLLISION  FREQUENCY  PROFILES 


Complete  electron  density  and  collision  frequency  profiles,  obtained 
from  the  analyses  described  in  Sections  IV-VII,  are  given  in  Figs.  46-53.  For 
the  purpose  of  establishing  error  bounds  on  the  plume  electrical  conductivity 
(a  function  of  n^  and  we  assume  that  the  collision  frequencies  have  been 
accurately  'computed  (via  the  collision  cross  sections  cited  in  Ref.  32)  and  all 
errors  are  "lumped"  into  the  error  bounds  placed  on  ng. 

Stage  1  ng  profiles  (Fig.  46)  are  estimated  to  be  correct  to  within  a 
factor  of  3,  up  to  the  position  of  the  first  Mach  disc  (xAj^  =  27).  Profiles 
downstream  of  this  point,  'within  the  second  shock  cell  are  based  on  the  center - 
lino  ng  profiles  (shown  in  Fig.  45),  which  have  been  arbitrarily  reduced  by 
about  an  order  of  magnitude  to  account  for  mixing  across  the  slip  line  (Section 
VI).  Values  of  n^  between  x/rj^  »  27  and.x  s  50  are  estimated  to  be  correct 
to  within  a  factor  of  10. 

The  Stage  2,  118  kft  plume  has  a  negligibly  small  Mach  disc  (Fig.  35), 
Thus,  conditions  downstream  of  the  shock  intersection  with  the  plume  center - 
line  can  be  determined  via  the  method  of  characteristics.  We  therefore  expect 
the  Stage  2  (118  kft)  values  of  n^  (Figs.  48a  and  48b)  to  be  accurate  to  within  a 
factor  of  3  for  x/v^  £  20,  and  within  a  factor  of  t  .out  6  for  x/rj^  >  20.  The 
dashed  curves  on  Figs.  48a  and  48b,  i.  c.  "slow  mixing"  were  determined  using 
the  original  Donaldson/Gray  eddy  viscosity  formulation.^  These  results  imply 
negligibly  small  mixing  rates  for  x/r^  <  *  30.*  Use  of  the  Tufts/Smoot  correc¬ 
tion’^  gives  more  realistic  turbulent  mixing  rates;  the  resulting  n^  profiles  are 
the  solid  curves  in  Figs,  48a  and  48b.  These  curvr:s  should  be  used  for  future 
electrical  conductivity  calculatic  ns. 

Stage  2,  200  kft  profiles  (Figs.  50  and  51)  have  been  computed  up  to 
the  approximate  position  of  the  first  Mach  disc.  The  dashed  curves  give  the 
results  for  "slow  mixing".*  Thus,  we  observe  that  at  x/rj^  =  60  the  inner  mixing 
boundary  is  located  at  r/x,^  *  12,  still  . far  from  the  plume  centerline.  These 
profiles  are  expected  to  be  accurate  to  within  a  factor  of  3,  up  to  x/r^  »  74. 
Calculations  downstream  of  x/rj^  s  74  must  include  an  analysis  of  flow  behind 
the  Mach  disc,  which  wm  made  in  the  present  study.  We  would, however, 
expect  increases  in  n^  siihilir  to  those  observed  in  the  Stage  I  plume  (Fig.  46). 

i 

rwrnnm  I . . .  y 

*  ThttSf  in  this  regihi^/th^  ’'slow  mixing"  psiofiles  are  identical  to  tbeinviscid 

profiles, 
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The  influence  of  mixing  on  the  and  profiles  for  Stage  Z  (300  kft) 
is  neglibibly  small.  Thus,  the  profiles  given  in  Figs.  52  and  53  are  those 
determined  from  the  inviscid  plume  calculations  only.  Values  of  n^  are 
expected  to  be  accurate  to  within  a  factor  of  3. 


IX.  SUMMARY  AND  RECOMMENDATIONS 

The  electron  density  and  collision  frequency  profiles  presented  here¬ 
in  have  been  determined  from  the  most  comprehensive  study  made  to  date  of 
Minuteman  exhaust  plume  electrical  properties.  Estimated  error  bounds  on 
computed  electron  densities  range  from  a  factor  of  3  to  10,  depending  on 
plume  location  and  flight  conditions.  Collision  frequencies  are  assumed  .to  be 
essentially  correct. 

In  the  process  of  computing  plume  electrical  properties  we  have 
shown  that: 

(1)  Gas  property  distrihutians  at  the  nozzle  exit  plane  depend  strongly  on 
the  weight  percent  61  particles  in  the  floW'-but  are  weakly  dependent 
on  the  precise  particle  distribution  ^within  the  range  of  2-8  p 
particle  diameters). 

(2)  Variation  of  the  important  rate  coefficients  for  the  charged  species 
reactions  within  their  ranges  of  uncertainty  do  not  have  a  significant 
influence  on  nozzle  exit  plane  electron  densities. 

(3)  Inviscid  plume  electron  densities  are  very  sensitive  to  the  rate 
coefficient  for  H  +  Cr~*’  HCl  +  e"  ,  which  produces  electrons  at 
relatively  low  temperatures. 

(4)  Electron  production  behind  the  first  Mach  disc,  and  within  the  second 
shock  cell,  is  significantly  high  for  Stage  1  (at  100  kft)  and  must  be 
accounted  for  in  determining  electron  density  profiles. 

We  reconunend  that  future  calculations  on  Minuteman  exhaust  plumes 
include:  (1)  the  influence.. of  heterogeneous  electron/positivo  ion  recombination 
on  plume  electron  densities,  (u)  the  use  coupled  ixMthod  of  characteristics/ 
turbulent  mixiag/nonequilibcium  chemistry  atU^lys is  for  those  flight  conditions 
where  the  n^xtog  region /is  si|nificahtly  luge  compared  to  the  plume  boundary 
radius  and  (iii)  ii^lS.r^ntstipjq  of  UcWique.  iMtiated  herein  to  determine 
flow  prqyerties  eptfj^si^idn  b^lj^W  disc  via  an  analysis  of 

'mixin^/noi^qqilibriumictisnUs^ry  ii|lohg;^a  slip  lii^ 
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TABLE  1 

•  f 

MINUTEMAW  PROPELLANT  COMPOBITIONS 
AND  NOZZLE  PARAMETERS  (U) 


Compoiition,  wt.  % 

Ammonium  Perchlorate 

Aluminum 

Binder 


8tage>  1 

(Thiokol  TPH-1011) 


Stage  2 

(Aerojet  ANB-30^^1 


70  73 

16  ,15 

14  (PBD/acrylic  acid  15  (PBCT) 

acrylinitrile) 


Elemente,  gram-atonia 

C 

H 

P 

N 

A1 

Cl 


0.929 

3.68 

2.47 

0.627 

0.593 

0.600 


0.865 

3.88 

2.50 

0.625 

0.556 

0.621 


Impuritiet,  ppm 

Na 

K 


180  180 

50  50 


Heat  of  Formation  (PropeUant), 
kcal/ mole 

-43.6 

Chamber  Pressure,  psia 

741 

Throat  Diameter,  in. 

1,ZZ 

Exit  Diameter,  in. 

22.8 

10 

rj^,  in. 

--- 

^'N)eq. , 

45.6 

•43.8 

445 

9.6 

48 

25 

48 


iJr 
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TABLE  IV 

ONE  DIMENSIONAL  EQUILIBRIUM  NOZZLE  EXIT 
CONDITIONS 


Stage  1  Stage  2 


/ 

Pressure,  psia 

11 

2. 1 

Temperature,  K 

2100 

1800 

Velocity,  ft/sec 

8500 

9300 

Mach  Number* 

2.86 

3.40 

a 

:v 


1.19,  W  =  28  for  equilibrium  gas/particle  mixture 
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TABLE  VI 

INPUT  DATA  FOR  2-PHASE  MOG  NOZZLE  CALCULATIONS 


« 

Stage  1 

Stage  2 

Chamber  Pressure,  psia 

741 

445 

Chamber  Temperature,  K 

3410 

3410 

Gas  Properties  (assumed  constant) 

1.27 

1.2'i 

^"g 

20 

20.  8 

Cp,  ,  Cc.l/mole-  K 
-  g 

Pr 

9.35 

9.5C 

0.78 

0.7J 

Ibrn/lt-sec 

6.0  X  10'5 

6.0  X  10-= 

n,  (p.  T'^) 

0.67 

0.67 

Particle  Properties  (AI2O3) 

Vp/Yg 

Diameter  (ii.)/weight  %  of  total 

0.393 

0.36 

2/10 

2/lC 

4/40 

4/4C 

6/40 

6/4C 

8/10 

8/lC 

Solidification  Temperature,  K 

2320 

2320 

Density,  Ibm/ft^ 

Specific  Heat  (T  =  2320  K) 

250 

250 

Cn  ,  cal/mole-  K 

34.6 

34.6 

C„  ,  cal/mole-  K 

Enthalpy  (  T  =  2320  K) 

33.0 

33.0 

h^,  kcal/mole 

87.8 

87.8 

h  ,  kcal/mole 

0 

61.8 

61.8 

Nozzle  Properties 

Inlet  Half  Angle,  deg 

40 

40 

Throat  Radius  of  Curvature,  in. 

7.99 

9.64 

*  at  combuatioii  chamber  temperature 
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TABLE  Vn 

INITIAL  CONDITIONS  FOR  NONEQUILIBRIUM 
NOZZLE  STREAMLINE  CALCULATIONS 


Stage  1 _  Stage  2 


Axial  Distance  from 


Throat,  in. 

-8.0^(C.L.) 

-8.4 

-7.0-  (wall) 

-5.9 

Pressure,  psia 

710 

400 

Temperature,  K 

3400 

3400 

Velocity,  fps 

900 

1070 

Composition,  mole  fraction 

CO 

2.57{-l) 

2.20(-l) 

CO2 

1.69(-2) 

2.22(-2) 

HCl 

1.45{-1) 

1.46(-1) 

H 

3.50(-2) 

4.20(-2) 

Ha 

2.93(-l) 

2.44(-l) 

OH 

8.20(-2) 

1.42(-2) 

HaO 

1.52(-1) 

1.93(-1) 

Na 

9.06(-2) 

8.70(-2) 

0 

6.30(-4) 

1.50(-3) 

Oa 

1.30(-4) 

4.50(-3) 

K 

1.22(-6) 

1.21(-6) 

KCl 

2.78(-5) 

2.74(-5) 

Na 

2.00(-5) 

2.15(-5) 

'  NaCl 

l.90(-4) 

1.92(-4) 

Cl 

1.20(-2) 

1.65(-2) 

Na+ 

8.50(-6) 

9.70(-6) 

K+ 

7.80(-6) 

8.70(-6) 

Cl" 

1.60(-5) 

1.80(-5) 

e" 

2.90(-7) 

4.00(-7) 

a 


negative  distances  are  measured  upstream  of  throat 
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TABLE  Vm 

REACTION  MECHANISM  AND  RATE  COEFFICIENTS^ 

FOR  NONEQUILIBRIUM  NOZZLE  AND  PLUME  CALCULATIONS 

k  =  AT'^exp  B/RT 


Error  Bound 


Reaction 

Ab 

N 

•  B 

Upper 

Lower 

I. 

O  +  O  +  M 

Oz  +  M 

l(-29) 

1 

0 

30 

30 

2. 

O  +  H  +  M 

OH  +  M 

l(-29) 

1 

0 

30 

30 

3. 

H  +  H  +  M 

Hz  +  M 

5(-29) 

1 

0 

30 

30 

4. 

H  +  OH  +  M 

— 

HzO  +  M 

2(-28) 

1 

0 

10 

10 

5. 

CO  +  0  +  M 

— 

COz  +  M 

l(-29) 

1 

-2,500 

3 

30 

6. 

OH  +  OH 

-* 

HzO  +  O 

l.O(-ll) 

0 

-780 

5 

5 

7. 

OH  +  Hz 

-*■ 

HzO  +  H 

3.6(-ll) 

0 

-5. 200 

3 

3 

8. 

O  +  Hz 

-* 

OH  +  H 

2.9(-ll) 

0 

-9,460 

3 

3 

9. 

H  +  Oz 

— 

OH  +  O 

3.7(-10) 

0 

-16,800 

3 

3 

10. 

CO  +  OH 

COz  +  H 

9(-13) 

0 

-1,080 

5 

5 

11. 

H  +  Cl  +  M 

— 

HCI  +  M 

4(.26) 

2 

0 

100 

100 

12. 

Cl  +  Hz 

— 

HCI  +H 

4(-ll) 

0 

-4,400 

10 

10 

13. 

HzO  +  Cl 

-* 

HCI  +  OH 

5(-ll) 

0 

-19,000 

30 

30 

14. 

OH  +  Cl 

— 

HCI  +  0 

3(-ll) 

0 

-5,000 

30 

30 

15. 

K  +  HCi 

KCl  +  H 

6(-10) 

0 

-5,000 

10 

30 

16. 

Na  +  HCI 

NaCl  +  H 

5(-10) 

0 

-8, 000 

10 

30 

17, 

K+  +  e'  +  M 

K  +  M 

2(-22) 

1.5 

0 

5 

5 

18. 

Na’'"  +  e"  +  M 

— 

Na  +  M 

i:5(-20) 

2 

0 

5 

5 

19. 

+  cr 

-• 

K  +  Cl 

l(-8) 

0.5 

0 

10 

100 

20. 

Na'*'  +  Cr 

Na  +  Cl 

3(-8) 

0.5 

0 

10 

10 

21. 

Cl  +  e  +  M 

-*• 

cr  +  M 

3(-30) 

0 

0 

30 

30 

22.^^ 

HCI  +  e- 

H  +  Cl’ 

1(.8) 

0 

-20,000 

30 

100 

Taken  from  Refs.  16  and  17,  whic.1  fully  discuss  the  sources  of  all  rate  coefficient 
data. 

^  cm-molecule-sec  units 

^  Estimated  values  (taken  from  Refs.  16  and  17  ,  which  give  precise  definitions  of 
error  bounds),  A  is  multiplied  by,  the  upper  bound  to  give  the  maximum  value 
of  of  k;  A  is  divided  by  the  lower  bound  to  give  the  minimum  value  of  k. 

Lower  bound  for  kgx  used  in  "standard"  plume  streamline  and  mixing/chemistry 
calculations  (see  Section  Vi)> 


TABLE  IX 


INPUT  DATA 
FOR  EXHAUST  PLUME 
METHODS  OF  CHARACTERISTICS  PROGRAM 


Altitude,  k£t 

Free  Stream  Mach  Number 

Stagnation  Pressure,  psia 

Stagnation  \^emp6rature,  K 

Specific  Heat  Ratio 

Molecular  Weight 

Mach  Number^’  ^ 

Centerline 
Nozzle  Lip 

Exit  Plane  Flow  Angle 
Nozzle  Lip,  deg 


Stage  1 


3440 


3380 

1.28 

21 


3.  22 
3.  22 


16 


Initial  data  line  is  nozzle  exit  plane.  Mach  number  and 
flow  angle  distributions  are  give  in  Fig.  18  (Stage  1). 

Con’':al  source  flow  in  the  nozzle  was  assumed  to 
establish  the  initial  data  line  for  Stage  2. 
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TABLE  X 

INPUT  DATA  FOR  MlXlNG/NONEQUILpRlUM 
CHEMISTRY  CALCULATIONS* 


Altitude,  kft 

Stage  1 

Stage  2 

100 

118 

200 

Jet  Velocity,  fps 

10000 

9300 

10150 

l- 

Jet  Temperature, °  K 

1200 

1170 

735 

Initial  Boundary  Pressure,  atm 

0,114 

0. 

0453 

0.0052 

Free  Stream  Velocity,  fps 

4850 

5300 

3850 

Free  Stream  Temperature,  K 

430 

430 

1400 

Pressure  distributions  are  given  in  Fig.  41, 

Jet  velocity  and  temperature  refer  to  rocket  exhaust  conditions  after  a 
Prandtl  Meyer  expansion  to  initial  boundary  pressure. 
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TABLE  XI 

INITIAL  CONDITIONS  FOR  NONEQOILIBRIUM  CHEMISTRY 
CALCULATIONS  DOWNSTREAM  OF  MACH  DISC 


Temperature,  K 
Pressure,  psia 
Velocity,  fps 

Composition,  Mole  Fraction 

CO 

CO* 

HCl 

H 

H* 

OH 

H*0 

Na 

O 

O* 

K 

KCl 

Nri 

NaCl 

Cl 

.Na+ 

K+ 

cr 


stage  1,  100  kft _ 

3330  (1060  K-free  stream)^ 
0.24 

1500  (9890  fps -free  stream)®^ 


2.35(-l) 
3.68(-2) 
1.58(-1) 
1.89(-3) 
2.82(-l) 
1.00 (-4) 
2.16(-1) 
9.10(-2) 
4.95(-7) 
2.00(-7) 
1.28(-7) 
3.55(-5) 
2.56(--6) 
2.20(-4) 
8.40(-4). 
1.87(-7) 
1.28(-6) 
L46(.6) 
4.34(.9) 


Refers  to  free  stream  tempefa^re  and  velocity  input  to  mix- 
ing/chemiitry  program.  %ecies  mole  fractions  in  free  stream 
were  assumed  equal  to  initial  species  mole  fractions  immediately 
downstream  of vMach  dU  c . 
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72-56 


Initi&l  Une  ««  for  Method  of  Ch«r»cteristic<  (MOC)  c^lculetion 


FIG.  6  GAS  TEMPERATURE  DISTRIBUTIONS  IN  STAGE  2  NOZZLE 


PARTICLE  VELOCITY,  fp 


4© 


FIG.  10  PARTIGLS  VELOCITY  DISTRIBUTIONS  IN  STAGE  2  NOZZLE 


PARTICLE  DIAMETER 


PARTICLE  TEMPERATURE 


FIG.  12  PARTICLE  TEMPERATURE  DISTRIBUTIONS  IN  STAGE  2  NOZZLE 


PARTICLE  DIAM 


GAS  VELOCITY,  fps 


FIG.  17  INFLUENCE  OF  PARTICLE  SIZE  ON  CENTERLINE  GAS 
.  VELOCITIES  IN  STAGE  2  NOZ Z LE 
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FLOW  ANGLE,  radians 


72-93 


DISTANCE  FRO^  CENTERilNE,  f/r^ 

FIG.  16  STAGE  I  mZZhE^mt  GAS  PROPERTIES 

Flow  angle  i  eagle  '\>p^pptSi  veloctty  vocler  eja4^^ino  axis 


r3.2 


GC 

UJ 

2.8  I 

3 

Z 

X 

•24  ^ 

S 


2.0 


,  >'  r) 


PARTICLE  TEMPERATURE 


72- »7 


rXG.  IH  STAGE  1  KOEZLB  EXlt  PLANE  PARTICLE  PROPERTIES 
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DISTANCE  FROM  CENTERLINE,  r/r^ 


FIG.  1 9b  STAGE  1  NOZZI.E  EXIT  PI^NB  PARTICLE  PROPERTIES 
Flow  amle  ?  »»gle  between  vi^tocity  vocto?  »,tiA  plume  axle 


FLOW  ANGLE.  RADIANS 


Fia.  20  SXAGB.2.  nozzle  ?;XIT  PLANE  GAS  PROPERTIES 

Flow  -  angle  between  velocity  vector  and  plume  axie 
1-0  Eqnllibrlum  P,  T,  M  from  Table  IV 
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MACH  NUMBER  PRESSURE,  psia 


DISTANCE  FROM  CENTERLINE  ,r/r^ 


FIG.  21*  STAGE  2  NOZZLE  EXIT  PLANE  PARTICLE  PROPERTIES 


0  0.2  0.4  0.6  0.6  t.O  1.2  1.4 

DISTANCE  FROM  CENTERLINE,  r/r,^ 


FIG.  2 lb  STAGE  2  NOZZLE  EXIT  PLANE  PARTICLE  PROPERTIES 
Flow  angle  =  angle  between  velocity  vector  and  plume  axis 
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PRESSURE,  psia  GAS  TEMPERATURE 


FIG.  22  INFLUENCE  OF  PARTICLE  SIZE  ON  STAGE  2 
NOZZLE  EXIT  PLANE  GAS  PROPERTIES 
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FIG.  24  H,  OH,  O,  Oj  MOLE  FRACTION  DISTRIBUTIONS 

IN  STAGE  1  NOZZLE 
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ELECTRON  MOLE  FRACTION 


AXIAL  DISTANCE  FROM  THROAT, inches 


FIG.  31  INFLUENCE  OF  RATE  COEFFICIENTS  ON  CENTERLINE 
ELECTRON  MOLE  FRACTIONS  IN  STAGE  I  NOZZLE 
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FIG.  33  PLUME  BOUNDARIES,  SHOCK  STRUCTURE 
AND  CONSTANT  MACH  NUMBER  CONTOURS 
FOR  STAGE  1  PLUME  (100  kft) 


woud  aoNvisia 
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FIG.  35  PLUME  BOUNDARIES,  SHOCK  STRUCTURE 
AND  CONSTANT  MACH  .NUMBER  CONTOURS 
FOR  STAGE  I  PLUME  (,118  kft) 


STAGE  2 
200  KFT 


.  36  PLUME  BOUNDARIES.  SHOCK  STRUCTURE 
AND  CONSTANT  MACH  NUMBER  CONTOURS 
FOR  STAGE  2  PLUME  (200  kit) 
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t  C  i  KON  DLNSITY 


ELECTRON  MOLE  FRACTION 


FIG.  39  ELECTRON  MOLE  FRACTION  DISTRIBUTIONS 
'  IN  STAGE  2  PLUME 
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ELECTRON  DENSITY,  m! 


PRESSURE,  psia 


FIG.  41  PRESSURE  DISTRIBUTIONS  ALONG  PLUME  BOUNDARIES 
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DISTANCE  FROM  CENTERLINE,  r/r 


DISTANCE  FROM  EXIT  PLANE  ,  K/r^ 


FIG.  42  DIVIDING  STREAMLINE  LOCATIONS 
Computed  from  Mixing  with  Nonequilibrium  Chemistry  program 


ELECTRON  MOLE  FRACTION 


i'Ui.  U  ITAIPERA'ItJRE  AND  ELECTRON  MOLE  EKAC  I  ION 

DIS  nUIiU  rtONS  along  dividing  si  ueamlines 

kzj/lOO 
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TEMPERATURE 


PRESSURE,  PS! 


DISTANCE  DOWNSTREAM  OF  MACH  DISC, 


FUi.  H  FS  I  IMA'I  ED  l-'KOl-’ER  T  Y  DIS'I  KlHU'l  IONS  DOW  NS  I  1(  EAM 
OF  MACH  DISC  -  STAGE  I  FEUME 


VELOCITY,  fps 


ELECTRON  MOLE  FRACTION 


FICi.  45  ELECTRON  MOLE  FRACTION  DISTRIBUTIONS 
DOWNSTREAM  OF  MACH  DISC  -  STAGE  I  PLUME 
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ELECTRON  DENSITY,  ml’ 
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DISTANCE  FROM  CENTERLINE. 


DISTANCE  FROM  EXIT  PLANE,  x/ ^ 
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ELECTRON  DENSITY,  ml 


I-'IG.  4aa  ELECTRON  DENSHY  PROFILES 
IN  STAGE  I  PLUME  (1 18  kft) 

0  <  x/rjvj<  45 
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FIG.  4Bb  El.ECTROy  DEXSITY  PROFIEES 
IN  Si  ACE  2  PEU.ME  lib  kit 


DISTANCE  FROM  EXIT  PLANE, 


FIG.  52  ELECTRON  DENSITY  PROFILES  IN  STAGE  2  PLUME  (300  kft) 
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APPENDIX  A 

NOZZLE  FLOW  CALCULATIONS 


Possible  approaches  to  computing  flow  properties  and  composition  in 
the  Minuteman  first  and  second  stage  nozzles  and  their  limitations  are  listed 
below  in  order  of  increasing  accuracy. 

(1)  One -Dimensional  Gas/Particle  and  Thermochemical  Equilibrium 
Neglecting  both  particle  thermal  and  velocity  lags  and  non¬ 
equilibrium  chemistry  effects  during  the  nozzle  expansion  makes  the 
subsequent  calculations  of  plume  electron  densities  meaningless. 

(2)  One -Dimensional,  Two-Phase  Flow  -  Equilibrium  or  Frozen 
Chemistry  -  Neglecting  the  2-dimensional  nature  of  the  flow  is  likely 
to  result  in  serious  errors  in  nozzle  exit  plume  properties.  Assum¬ 
ing  frozen  electron  mole  fractions  in  the  nozzle  will  drastically  over¬ 
estimate  exit  plane  electron  densities;  assuming  equilibrium  flow  will 
drastically  underestimate  exit  plane  electron  densities. 

(3)  Two-Dimensional,  Two-Phase  Flow  -  Nonequilibrium  Chemistry 
Along  Streamlines  -  In  this  approach  the  two  phase  nature  of  the  flow 
is  accounted  for  together  with  the  two-dimensional  effects.  The  influ¬ 
ence  of  nonequilibrium  chemistry  on  the  streamline  locations,  pres¬ 
sure,  and  gas  and  particle  temperature  and  velocity  is  neglected  (a 
reasonably  good  assumption).  The  above  properties  are  input  to  a 
nonequilibrium  calculation  of  neutral  and  charged  species  concentra¬ 
tions  along  streamlines. 

(4)  Two  -Dimensional,  Two-Phase.  Flow  Coupled  with  Nonequilibrium 
Chemistry  -  A  program  which  incorporates  the  influence  of  non¬ 
equilibrium  chemistry  on  pressure  and  gas  and  particle  temperature, 
and  velocity  yields  an  "exact"  solution  to  this  problem. 

After  consulting  wirh  those  research  laboratories*  likely  to  have 
program  (4)  we  concluded  that  none  was  available  for  use  in  this  study.  We 
therefore  chose  analysis  (3)  to  make  the  Minuteman  nozzle  calculations. 


Including  United  Aircraft,  Lockheed -Huntsville,  Grumman,  Purdue 
University 
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Results  of  a  comprehensive  investigation  of  Minuteman 
exhaust  plume  electrical  properties  are  presented.  Electron 
density  and  electron -neutral  collision  frequency  profiles  are  given 
for  the  Stage  1  plume  at  100  kft  and  the  Stage  2  plume  at  118,200 
and  300  kft.  This  work  is  motivated  by  the  need  to  determine  plume 
electrical  conductivities,  input  data  to  an  evaluation  of  EMP/exhaust 
plume  interactions.  The  calculations  are  described  in  detail  and 
include  analyses  of:  (i)  two-phase  nozzle  flow  via  the  method  of 
characteristics,  (ii)  nonequilibrium  chemistry  along  nozzle  and 
plume  streamlines,  (iii)  inviscid  plume  properties  via  the  method  of 
characteristics,  (iv)  turbulent  mixing  /  nonequilibrium  chemistry 
along  curved  plume  boundaries,  and  (v)  nonequilibrium  chemistry 
behind  the  Mach  disc.  The  results  include:  gas  and  particle  prop¬ 
erties  and  cheznical  kinetic  distributions  of  all  neutral  and  charged 
species  concentrations  in  the  nozzle:  inviscid  boundaries,  shock 
structure,  Mach  number  contours  and  electron  densities  in  the 
pldtne* 


.  .  Rui  .  .  .  . 

mm 


JlMCiiMSlFXSD 


